Numerically efficient version of the T-matrix method. 
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Abstract 

CN ' A version of the projection method for solving the scattering problem for acous- 

tic and electromagnetic waves is proposed and shown to be more efficient numeri- 
. cally than the earlier ones. 

O ' 
m . 

§ : 1 Introduction 

o : 

Consider the scattering problem: 

(V 2 + A; 2 ) u = in D' := IR 2 \D, (1.1) 

J* i u = on S := 3D, (1.2) 

u = e **«-* + Vj (1.3) 

where D is a bounded domain in IR 2 with a piecewise smooth boundary S, and the 
' scattered field v has the following asymptotics: 

e ir ( 1 \ x 
ctf ■ u = A(a',a,fc)— = + o -= , r = Id -»• oo, - = a' <E S 1 , (1.4) 

r \ vT / r 



and S" 1 is the unit sphere in M 2 . The coefficient A is called the scattering amplitude, 
k = const > is fixed, and u := u(x, a, k) is called the scattering solution. 

Let g(x, y) := jHq(Jc\x — y\). Take k — 1 in what follows without loss of generality. 
Then 



e .-iy°-x i _ i 1 \ i.. i _ _ ..o . y 



(V 2 + l)g = -S(x - y) in R 2 , g = c^=e~ iv ' x + o — = as Id -> oo, y 

vlil VvW ' 12/ 



(1.5) 



const = ttt^- One has, using Green's formula: 



/ g(x, s)un(s' ,a)ds' , x G Z?' , fe = 1, (1.6) 



where is the normal derivative of u on S, N is the exterior unit normal to S. 
Taking x E D' to S and using (1.2), one gets 



Th:=f, h:=u N , f:=e ia - s , k = 1, (1.7) 
T/i:= / g(s,s')h{s')ds'. (1.8) 



JS 

Equation (1.7) is the basic equation studied in this paper. 

The T-matrix method is described in [0], and analyzed in jj]] mathematically. 

The purpose of this paper is to give a version of this method for solving equation 
(1.7), and to analyze this version from the computational points of view. This is done in 
sections 2 and 3. 

Let H l = H £ (S) be the Sobolev spaces, I = 0, 1, and 

Qh:= f -L\ n —h(s')ds', (1.9) 
Js 2tt r ss , 

where a = const = diam D > 0, so that inf SjS / g s > 1. 
Write (1.7) as 



Qh + Kh = f,fe H\ Kh := -— / h(s')ds'+ [ 

2?r Js Js 



27r r.«/ 



/idy (1.10) 



The operator A = Q~ X K is compact in H 1 , and the operator Q is an isomorphism 
between H° and H 1 (see Lemmas 1.1 and 1.2 below). Moreover, Q is a selfadjoint 
compact positive operator, in H° : (Qu,u) > if it 7^ 0. In Remark 2.1 below we show 
a possible usage of (1.10). Let {<Pj} 3=1,2,... be a Riesz basis of H°, that is, every element 
u G H° is uniquely representable as a convergent in if series 



and 



00 00 
raV] I c j 1 2 < ||u||q < mV |cj| 2 ,0 < m < M, m,M= const > 0. (1.12) 
j=i i=i 

Let us prove the following Lemma 1.1: 

Lemma 1.1. If k 2 is not a Dirichlet eigenvalue of the Laplacian in D then the operator 
T defined by (1.8) is an isomorphism of H° onto H 1 . 
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Proof. This result is established in [JTJ] so we only indicate the basic points of the proof. 
The assumption of Lemma 2.1 implies the injectivity of T : if Th = then the function 
w(x) := J s g(x, s')h(s')ds' solves the Dirichlet problem for the Helmholtz operator in 
D and in D', and satisfies the radiation condition at infinity. Thus w — in D by 
the assumption of Lemma 2.1, and w = in D' by lemma 1 in Q, p. 25. Therefore 
h = u~n — ujf = u, where we have used the jump relation for the normal derivative of the 
single-layer potentials. The operator T : H° — > H l is of Fredholm-type: it can be written 
as T = Q+, K where Q : H° — > H 1 is an isomorphism and K is compact as an operator 
from H° into H 1 . The injectivity of T together with its Fredholm property imply the 
conclusion of Lemma 1.1. □ 

We assume throughout this paper that k 2 = 1 is not a Dirichlet eigenvalue of the 
Laplacian in D. 

In the literature @ one usually means by the T-matrix approach (in acoustic and 
electromagnetic wave scattering theory) a projection method for solving equations of the 
type (1.7) with the following choices of the basis functions: ip m = e %mb ' H m (kr(6)) or 
ip m = e me J m (kr(6)). 

These choices lead to the following difficulties discussed in jlj: as the number J of 
these functions grows: J — > oo, the condition number of matrix in (2.2) (see below) 
grows exponentially and depends strongly on the geometry of S. In contrast, in our 
version of the method, the condition number of a^- remains bounded as J — > oo. 

Lemma 1.2. The operator Q~ X K is compact in H 1 . 

Proof. The kernel of the operator K, defined by (1.10), and its first derivatives are 
continuous functions of s and s' running through bounded sets, including the diagonal 
s = s'. By Lemma 1.1 the action of Q _1 is equivalent (up to the terms preserving 
smoothness) to taking the first order derivatives. Therefore the conclusion of Lemma 2.1 
follows. □ 



Remark 1.1. Let us outline a possible way to use the splitting in equation (1.10). The 
idea is simple: Q > c > is positive definite. Write (1.10) as 

h + Ah = F, F:Q- X f. (1.13) 

Inverting Q numerically is a relatively easy problem since Q is positive definite. Using an 
orthonormal basis {(pj} of H°, one can write the projection method for (1.13), namely: 

j 

h J = 

i=i 

J 

hi + Y^ A ii h i = l<i<J, (!- 14 ) 

where := (Aipj,ipi). 
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The matrix in (1.14) is = 5^ + A^. If ||A|| < 1, then system (1.4) can be solved 
by iterations numerically efficiently This happens if ka << 1, where a = diam D, but 
it may happen when the above condition does not hold. 



2 Solution of the basic equation 

Let us look for an approximate solution to (1.7) 



j 

h = S Z2 c i¥j, (2-1) 

i=i 

where ipj is a Riesz basis of H°, and 

j 

a n c 3 f- l<i<J, (2-2) 

3=1 

where 

/i := (f,Tipi), Oij := (T^-,T^)o, (2.3) 

and T is defined in (1.8). Since T is injective, the elements {TV-,} are linearly indepen- 
dent. Therefore 

det(a ij )i<ij<j ^ VJ = 1,2,... (2.4) 



Definition 2.1. ^4 system {ipj} is a Riesz basis of a Hilbert space H if there is an 
isomorphism B of H onto H such that Bipj = ej, where {ej} is an orthonormal basis of 
H. 

One gets system (2.2) by solving the following minimization problem: 

j 

\\^2 c J Tl Pj - /111 = min. (2.5) 

3=1 

Denote by {cf ^}i<j<j the unique solution to (2.2) or, equivalently, to (2.5). 

Lemma 2.1. {Ttpj}i<j <00 is a Riesz basis of H 1 if the system {<Pj}i<j<<x is a Riesz basis 
ofH°. 

Proof. Let / G H 1 be an arbitrary element of H 1 . Denote T~ 1 f := h E H°. 
Since {<fij} is a basis of H°, one has 

oo 
3=1 
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(ej, e i ) 1 = 5. 



If c j^ L Pj = 0, then, applying the continuous operator T 1 , one gets X^jli c jfj = 0, so 
Cj = for all j since {<£j}i<j<oc is a basis of if . We have proved that {T(pj}i<j <00 is a 
basis of if 1 . 

Let us prove that if { l fj}i<j<oo is a Riesz basis of if then {fVj}i<;<oo is a Riesz 
basis of if 1 , that is, there exists an isomorphism B of if 1 onto if 1 such that ST 1 ^ = e,-, 

1 z = j 

Let {e,} be an orthonormal basis of if 1 . Define a linear operator F : if 1 — > ff 1 by 
the formula: 

oo oo 

F J2 c 3 e i := J2 c j T( Pj> 

3=1 3=1 

in particular, Fej = Ttpj. Let us prove that F is an isomorphism of if 1 onto ff 1 . If this 
is proved, then f? := F" 1 , and Lemma 2.1 is proved. 

Clearly F is linear, is defined on all of if 1 , and is continuous. Only the continuity of 
F needs a proof. 

Let u n > u in if 1 . Then u n = J2j*Li ^ e 3^ u = J2j*Li c j e j' Yl'jLi \ c j^ ~ c j\ 2 °> 
*K = E°°=i c^Tipj, Fu = 52°° =1 Cj T Vj . Thus: 

oo 

(n) 



|F Mn -F W ||? = ||^( c W-c,)r^||? 



oo oo 

r - c >,n? < ii t ii 2 ii - ^Villo < 

oo 

WTfM^cf - Cj \ 2 ^0 asn^oo, 
j=i 

where we have used the assumption that {<^j}i<j<oo a Riesz basis of if . 
Thus F is a linear continuous, defined on all of if 1 operator. 
Therefore F is bounded. 

Let us check that F is injective: if m 6 ff 1 , u = Ejli c i e j> an d Fu = then 
^jLityTipj = 0. Apply T" 1 and get Yl'jLi CjVj = 0- Thus, Cj = Vj, since {</?,} is a 
basis. The injectivity of f 1 is proved. 

To complete the proof one has to check that the range of F is the whole space if 1 . Let 
us do this. Take an arbitrary / G ff 1 and define h := T~ 1 f G if . Let /i = Ejli c i^Pji 
then 

oo oo 

Th = f = c i T ¥j = F Y1 c J e 3- 

3-=l 3=1 

Therefore F is an isomorphism of if 1 onto if 1 , and {Tipj} is a Riesz basis of if 1 , as 
claimed. Lemma 2.1 is proved. □ 
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Let us summarize the proposed method for solving the basic equation (1.7): 
Step 1. Choose a Riesz basis { i Pj}i<j« x in H° = L 2 (S). 
We discuss this choice below. 

Step 2. Calculate the matrix entries and the numbers fi, 1 < i,j < J, where J 
is an a priori chosen integer. 

Step 3. Solve linear system (2.2) numerically. 

The matrix in (2.2) has condition number that remains bounded when J grows, as 
follows from Lemma 2.1. 

Let us discuss the choice of the basis {tpj}. 

Assume that r = r{6) is the equation of S in the two-dimensional case that is, S is 
star-shaped. The element of the arc length of S is ds = ■\/r" 2 {9) + r 2 (9)d9 := a(9)d9. 
Let S 1 denote the unit sphere S 1 := {x : x G M 2 , \x\ = 1}. Choose 

1 cos(m6') sin(m^) 

^ c °( s ) = ~lK=fm^ ¥cm{s):=^== tp sm {s) = ^==, m = l,2, ... (2.6) 
■ s /27ra(9) ■ s /7ra{9) yjira{9) 

where s = (r(9), 9). 
Then 



and 



1 f i 

(p cm (s)ip sm '(s)ds — — / cos(m9) sm(m'9)d9 = 0, 

vr Jo 

I (Pcm^cm'ds b~mm'-> I ^Psm^Psm'ds 0~mm' j 

Js Js 

so that {(p m (s)} is not only a Riesz basis of H° = L 2 (S), but an orthonormal basis of 
H°. 

Similar construction holds in R 3 , where the normalized spherical harmomics Y im (a) 
are used in place of cos(m6') and sm(m9), a = (9, ip) is the unit vector in R 3 , a e S 2 , S 2 
is the unit sphere in R 3 . 
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